arXiv: 1505.04194v4 [cond-mat.str-el] 9 Mar 2016 
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We address the problem of Dirac fermions interacting with transversal optical phonons. A gap in 
the spectrum of fermions leads to the emergence of the Chern-Simons excitations in the spectrum of 
phonons. We study the effect of those excitations on observable quantities: the phonon dispersion, 
the phonon spectral density, and the Hall conductivity. 
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I. INTRODUCTION 


During the last decade there has been consider¬ 
able progress in experimental studies of the two- 
dimensional (2D) Dirac gas in various fields. This 
was initiated by the investigation of exfoliated and 
epitaxial graphene l-Qj followed later by experi¬ 
ments with silicene ^ and with a number of chem¬ 
ical compound^commonly referred to as topolog¬ 
ical insulators [aQ. Other realizations are optical 
hexagonal lattices filled with ultracold atoms [1- 
fl^ and photonic crystals [ail. Besides the lin¬ 
ear Dirac spectrum with nodes, particle-particle 
interactions can play an important role, depending 
on the specific system. For instance, the coupling 
to lattice vibrations (phonons in condensed matter 
realization [i|-[l, shaking in optical lattices j silp 
and photonic crystals with Dirac spectrum [l3lll4| l 
can lead to instabilities flail^. Such an instabil¬ 


ity was recently observed in graphene [I3- The 
complexity of the phonon modes can give rise to 
the plethora of other interesting effects. This is 
related to the fact that in-plane optical phonons 
act like effective gauge fields that couple to the 
Dirac particles. Moreover, it has been known for a 
long time that the coupling of massive Dirac par¬ 
ticles to gauge fields leads to Chern-Simons ex¬ 
citations |18l - [^ . It is crucial to note that the 
Chern-Simons term is created by the expansion 
of the fermion determinant with respect to a vec¬ 
tor field, where the latter is usually a gauge field. 
We will show in this paper that an optical in-plane 
phonon mode, which is a massive vector field, can 
also create such a term. In condensed matter sys¬ 
tems there is no ambiguity from the UV regulariza¬ 
tion because of a natural momentum cut-off from 
the underlying lattice. The main objective of this 
paper is to study how the Chern-Simons term af¬ 
fects the transport properties. In particular, we 
will focus on the quantum Hall effect and study 
the Hall conductivity when the Dirac fermions are 
gapped. 


II. THE MODEL 

Our microscopic lattice model describes the in¬ 
teraction between spinless fermions and monochro¬ 
matic transversal optical phonons on a two- 
dimensional hexagonal lattice. The dynamics of 
the free fermionic quasiparticles is given by the 
tight-binding Hamiltonian 

Ho = -tY^ (cld^, + 4'Cr). (1) 

(rr'} 

where c and d denote fermionic species identified 
with each sublattice of the hexagonal lattice and 
t ^ 2.8eV is the hopping amplitude. We use the 
unit system with lattice constant, elementary elec¬ 
tric charge, and H = 1. The index (rr') suggests 
the summation over next neighbors, while the ab¬ 
sence of the parenthesis means summation over 
each sublattice index separately. The Hamiltonian 
is readily diagonalized by the Fourier-transform 
and yields a well-known spectrum with two degen¬ 
erated parity symmetric Fermi points. As a strik¬ 
ing feature, close to those Fermi points the fermion 
dispersion is linear and therefore describes mass¬ 
less Dirac particles. In order to model the fermion 
mass we have to account either for a staggered po¬ 
tential altering its sign from one lattice site to an¬ 
other [l^, or for a spin-orbit coupling term 
The parity between the cones is broken by exposing 
the system to the chemical modification, e.g. by a 
periodic flux or a spin texture [1^. Combin¬ 
ing fermionic operators to spinors Tj. = (c^,d^)^ 
and = (c, d)J, the mass term reads 

Hl= ( 2 ) 

T r.r' 

where t' is the hopping amplitude between second 
nearest neighbors (t'/t ^ O.I 0 ) and the matrix 
element of the Haldane term is given by [ 2 ^ 

i=l,2,3 

with the adjustable Peierls phase 4>, and a^ denot¬ 
ing the positions of the second nearest neighbors 
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FIG. 1: (Color online) The fermion band structure of the considered model. The parity symmetry between both 
fermion copies is explicitly broken as visualized by the gap 2m ~ t at the Fermi point in the negative momentum 
-K versus the gapless state at the Fermi point at the positive momentum K= dvr/Sv^. The total band width 
2Ao ~ 6.5t is only slightly larger than that of the pristine tight-binding Hamiltonian in Eq. 0. 


on the hexagonal lattice. A particular choice of 
both parameters ^ = 3t'(l + v^)/v^ and (/> = 7r/4 
opens up a gap at one Dirac point with size 2m ~ 
2 X 0.75t, while keeping the other Dirac point gap¬ 
less. This band structure is visualized in Fig. [TJ 
Moreover, changing the value of the Peierls phase, 
it is also possible to get into the situation where all 
cones are gapped with different signs of the Dirac 
mass. For the chosen parameter set, all elements of 
the conductivity tensor of the Haldane model are 
nonzero [mi, where the gapped (gapless) chan¬ 
nel contributes to the Hall (longitudinal) conduc¬ 
tivity. But since we are primarily interested in 
the anomalous Hall conductivity, we focus on the 
gapped channel only. Recently, an experimental 
realization of the Haldane model was reported in 
Ref. [ 3 . 

In our model the monochromatic phonons ap¬ 
pear as fluctuating bonds between neighboring 
sites [2^-[3l|- The corresponding Holstein Hamil¬ 
tonian reads 

(rr') 


The operator 6 ,.^, (&^/^) destroys (creates) a 
phonon between sites r and r'. In the continuous 
limit the phonons couple to the spinors via com¬ 
binations of non-diagonal Pauli matrices ± i <^2 
with the coupling strength a. While the phonon 
frequency 7 ^ 160 — 170meV is intrinsic to the 
hexagonal lattice [l3111ll|, the electron-phonon 
coupling strength a, measured in units energy • 
length, can be varied in the experiment. 

In the continuous limit we use the low-energy 


approximation for massive Dirac fermions. This 
approximation is well known from the litera¬ 
ture [^1^ and we do not discuss it here. Next, 
we introduce the coherent state functional integral 
with the partition function given by 

Z = JV[ijj,ip,A] (5) 


with the 2+1-dimensional Euclidean action 


S[tlj,tlj,A] = -^ J d^x Al 

+ j d^x i)[$ + mA ( 6 ) 


where $ = + 7i9r,_+ > 7 = (7o, 7i, 72 ) = 

(ct3, -ZCTiCTa, -fcr2CT3), V’ = '0''’7o) ^nd / d^x = 
f drd^r, r denoting the imaginary time. The 
Fermi velocity vr = \/3t/2 is removed from the 
fermionic part by rescaling —>■ VpX^ (not r) and 
"0 —>■ '\l}lvp. Real valued components of the phonon 
field /i = 1, 2 are obtained via 




Al — f A2 Al + f A2 

V2a ’ 7^’ 


( 7 ) 


i.e., by changing from the U{1)- to the 0 ( 2 )- 
representation of the phonon sector. It is impor¬ 
tant to stress that, in contrast to gauge fields, our 
A^ are massive fluctuations of the lattice, which 
cannot be constrained by gauge fixing. Formally 
however, the classical equation of motion of the 
field A has approximately the form of the Lorenz 
gauge condition, cf. Appendix lAl 
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The coupling parameter g in Eq. m is related 
to the electron-phonon coupling a in ([3]) through 
2g = a‘^/{vpj). Finally, the slow lattice dynam¬ 
ics drAfj,) is neglected in comparison to the 
much faster electron dynamics. Integrating out 
the phonon fields yields an interaction term 
for the fermions {g/4:){'ij}j^ij})‘^. This resembles the 
standard Thirring current-current interaction di¬ 
ll^ with one crucial difference, though, that we 
have only two spatial currents here. Due to the 
anticommutativity of the Grassmann field, this 
interaction can be rewritten as the usual repul¬ 
sive Hubbard interaction {g/2){ip''^il;)'^. Through¬ 
out the subsequent calculations we implement a 
large but finite UV-cutoff regularization scheme 
with only spatial components of the frequency- 
momentum vector being regularized, while the in¬ 
tegration over the Matsubara frequencies stretches 
to infinity, hence only considering the T = 0 case 
here. Fig.[T]suggests the appropriate energy cutoff 
to be Ap - 3.2t - 10®K. 


analysis [3^ . Since the Dirac Hamiltonian is un¬ 
bounded, some terms in the perturbative expan¬ 
sion are plagued by ultraviolet divergences. There¬ 
fore, the fermionic determinant should be prop¬ 
erly regularized. There exist many regulariza¬ 
tion schemes diil, but in our case the spec¬ 
trum is naturally bounded by the band of width 
of 2 Aq ~ 6.5t. Therefore, all terms in the gra¬ 
dient expansion acquire small corrections; e.g., of 
linear or quadratic order in to/Aq for the gradient 
expansion. Higher order terms are negligible, as 
the analysis of the conductivity tensor for the full 
Haldane model in Ref. [2^ has shown. 

Performing a gradient expansion up to second 
order in the momenta, cf. Appendix [B1 we obtain 
for the effective action three terms as 

5[A] ^ -I- -I- Sj^^[A], (11) 

with the mass term 

5^ [A] =Ajd^xAl, (12) 


III. PHONON THEORY where 


After integrating the fermions in Eq. (j6|), we 
obtain an effective action in terms of the phonon 
fields A 


5[A] = 


^d^^-trlog 

2 g 


+ m + i 


A^ 

J 


( 8 ) 


where the functional trace tr consists of the three 
dimensional integration and the trace with respect 
to the Dirac space. Next we expand the tr log- 
term of Eq. (|5]) up to second order in the fields 


—tr log 


+ m+ ^A 7 


-tr log[^ -I- to] -^tr[^ -I- to] ^ a 7 

v 2 


--tr[^-bTO] ^A^7^[^-fTO] M^7^. 


(9) 


The linear term in this expansion is traceless, sig¬ 
naling that in classical approximation we get A^ = 
0. Thus, Gaussian fluctuations of the phonon field 
are relevant and lead to the effective action 

5[A] - ^ K 

-itr[^-bTO]"^A^7^[^-f to]"M^7^, (10) 


where higher order terms can be neglected be¬ 
cause they are irrelevant in terms of a scaling 


A 


1 _ ^0 
2g IOtt 


(13) 


and Aq denotes the UV-cutoff. The second term 
is the standard Chern-Simons term 


^CS i-^] 


se 




IOtt 


(PP 




(14) 


where denotes the antisymmetric tensor, s = 
sgn{m) and P = (P 01 P 11 P 2 )"'"; arising due to the 
energy gap in the spectrum of 2-1-1-dimensional 
Dirac fermions. We have neglected corrections 
C>(to^/Aq) in the prefactor of this term (cf. Ap¬ 
pendix [B]), assuming that the gap (fermion mass) 
is much smaller than the cutoff momentum [TsUToj] . 
The third term is a Maxwell-like term 




d^P ^ P^d^u - P^Pu ^ 

(27r)3 dSTrlTO] ’ 


( 15 ) 

although its structure is somewhat different from 
the usual Thirring model [Tsl - I^ . As a con¬ 
sequence of the missing gauge invariance of the 
phonon field, it cannot be written as a product 
of Maxwell field tensors = df^Aj^ — d^A^j^. The 
detailed derivation of Eqs. (fT^ - (|T5]) is given in Ap¬ 
pendix [BJ The presence of the term (flSl) in the ef¬ 
fective phonon action makes this action different 
from the otherwise similar Ghern-Simons-Proca 
action [ 3 ^ . 

Obviously, the phonon energy gap in Eq. 0 
can become negative, signaling an instability in the 
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functional integral. This instability is due to the 
transition into another phase [1^1^ . With the 
values Aq ^ 3.2t and m ^ 0.75t, which we have 
used for Fig.[l] this instability is irrelevant at t = 1 
for sufficiently small electron-phonon coupling: 

This inequality is valid in the low-energy approx¬ 
imation given by the action ([5|). 


to understand the role of the Chern-Simons term 
in the structure of elementary excitations we can 
formally set b = 0. In this limit we obtain 

e^ip) = Jp^ + -, e_ip) = J^, (23) 

V a V a 

which means that the e_-mode has a flat band and 
the energy gap between both modes disappears. 
Going even further and setting A = 0, leaves us 
with the linear mode of an acoustic phonon. 


IV. PHONON PROPAGATOR 


V. HALL CONDUCTIVITY 


The Gaussian action in Eqs. (dlHISl) enables 
us to calculate the correlations of the phonon 
fields without making any additional assump¬ 
tions. The calculation of the phonon propagator is 
straightforward and yields 

{A,{P)A^{-P)) = ^U^P), (17) 

where 

" a^[pl + el{p)][pl + el{p)] 

/ A -t a(po -f pf) sbpo + apiP2 \ 

X . (18) 

\ -sbpo + apiP2 A -h a(po +pl) ) 

Here b = I/IOtt denotes the Chern-Simons and 
a = l/(487r|TO|) the Maxwell weight factors. Large 
(small) values of \m\ favor the Chern-Simons 
(Maxwell-like) term. The elementary excitations 
of phonons are given by two gapped modes shown 
in Fig. [21 where the dispersion relations read 


e±(p) 

= Vfip)Pt(p)^ 

(19) 

. , 1 / 

Y A 62 \ 

(20) 

fip) = TT 

P^ + 2-1—^ , 

2 \ 

^ a J 

<(p) = A 

/(^-5)^45. 

(21) 


A physical quantity, which can be directly evalu¬ 
ated and also experimentally measured, is the Hall 
conductivity. It is calculated from the Kubo for¬ 
mula (^ : 

Oqr r 

= Y J e"“^((V’7^V^)„(?/’7^?/')o): (24) 


for fi V, where (• • ■) denotes the normalized func¬ 
tional integral with the action of Eq. m- Then in 
terms of the phonon field the correlation function 
in Eq. (HI becomes: 


- (25) 

Some details of this mapping are shown in Ap¬ 
pendix |Cj With the help of Eq. (fTSl) we obtain the 
Hall conductivity as 


2nb 

“ aV [w2 + e^(0)][u;2 + e^(0)]' 


(26) 


Being mainly interested in the behavior of the DC 
Hall plateaux, we send w —>■ 0. Since the pa¬ 
rameter s = sgn(m) allows values ±1, there are 
two plateaux at positive and negative mass value 
for each combination of /i and v. In the non¬ 
interacting case the distance between the plateaux 
is 2m. From e^(0)e?_(0) = jc? we obtain for 
the system with electron-phonon interaction 


^dc _ 

“ 8^2A2 ’ 


(27) 


and p2 = p 2 _|_ While for p <C 2A both 
modes grow ^ p^, at p » 2 A the features of non¬ 
interacting model are recovered (cf. Fig. |2|): the 
lower mode approaches the finite value 

[a W 

e_(p^oo) = y'- + ^, (22) 

while the upper mode crosses over into the linear 
regime of Dirac electrons e+(p ~ oo) ^ p. In order 


where the definition of b has been used. For weak 
interaction (or large mass) we approximate A ^ 
l/2g and reproduce the Hall conductivity of non¬ 
interacting massive Dirac fermion gas 

s 

For a finite band width we have p^A^ ^ 1/4, such 
that the dc plateaux are shifted away from each 
other by the electron-phonon interaction. 
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FIG. 2: (Color online) Left box: Phonon modes from Eq. (I19II normalized by e_(oo) depicted versus the momentum 
in units of the phonon mass, see the discussion in the main text. Right box: The phonon spectral density Eq. (p9l) 
calculated for the same parameter set. 


VI. PHONON SPECTRAL DENSITY 


The spectral density 

= -Im J fw + 77) (29) 


is experimentally accessible in neutron or Bragg 
scattering experiments. Below we discuss the di¬ 
agonal component The excitation modes 

shown in Fig. give a plausible hint to the ex¬ 
pected shape of the spectral function; It must com¬ 
prise two separate structures, corresponding to the 
two bands and the separation by the gap of size 
e^(0) — e_(oo). Moreover, it must reveal an addi¬ 
tional energy gap up to e_(0). By retaining only 
terms which are resonating at positive frequencies, 
we obtain 


od"”/ (2,)3X 

V + ^ . I 

u — — irj uj — e_ — IT]) 


(30) 


where the coefficients of the decomposition read 


A — a(e^ — 0.5p^) 

2e+(e^-e2_) 

A — a(e^ — 0.5p^) 

2e-(4-e-) 


(31) 

(32) 


Below we evaluate the spectral function for two 
extreme parametric regimes, corresponding to the 
very small and very large fermionic mass m. 


Chern-Simons regime: For large fermionic mass 
the Chern-Simons contribution is dominant. Then 
for very small a we may approximate A and B as 


lim = 

a ->-0 


1 


lim ^ = 0 . 

a -^0 


(33) 


After an integration over we obtain 


oCS 



Stto 



(34) 


Since h/a = 3|m|, the effect of the Chern-Simons 
term can be detected as an excitation continuum 
with frequencies larger than 3 |to|. Remarkably, 
this contribution is indifferent to the value of the 
phonon gap A. For comparison we show in Fig. |3] 
the full numerical evaluation of S^,^{ui) in the para¬ 
metric regime close to Eq. dSl- Actually, the ex¬ 
citation continuum in Eq. (IMl) is cut off at a scale 
related to the band width Aq. This is not included 
here. 

Maxwell-like regime: The contribution from 
the Maxwell-like term is dominant for a small 
fermionic mass. Then 


lim — = 


4a 


A 

a 


V ^ 1 

lim — = — 

b^o a 4\/aA 


(35) 


and the energy modes reduce to the expressions in 
Eq. (^5)) . Due to the flat band e_, the integration 
requires some care. In particular, one needs to 
keep r] nonzero, which ultimately leads to 

^ ^0 _ H _ 

^ 8v^(2^)2(cc-yA7^)^+^2 
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FIG. 3: (Color online) Left box: The spectral function Eq. (1291) close to the Chern-Simons regime. Blue (dashed) 
peak on the left shows the remnants of the resonance at e_ which disappears as o —^ 0. It is positioned at 
u! ~ A/a. Right box: The spectral function close to the Maxwell-like regime. The flat band and continuum 
contributions are shown in blue (dotted line) and red (solid line). The black (dashed) curve is the superposition 
of the both contributions. 


where the first sharply peaked contribution is a 
consequence of the flat band, cf. Fig. |31 Again, 
the 0-function which cuts off the excitation con¬ 
tinuum at the scale ~ Ag is not shown in Eq. dMl). 


VII. DISCUSSION AND CONCLUSIONS 

Our analysis of massive Dirac particles, which 
are coupled to optical in-plane phonons, has re¬ 
vealed an effective phonon model with a Chern- 
Simons term in Eq. (|T^ and a Maxwell-like term 
in Eq. m- These two terms compete, where the 
former is favored for a large Dirac mass. The op¬ 
tical phonon dispersion, which is flat without cou¬ 
pling to the Dirac particles, has two branches and 
becomes parabolic for small momenta (cf. Fig. 
[21). For large momenta one branch is linear and 
the other is flat. The spectral density of the 
two modes is shown for large Dirac mass (Chern- 
Simons regime) in Fig. [3] and for small Dirac mass 
(Maxwell-like regime) in Fig. |21 These results 
indicate a separation of the two modes which is 
obviously caused by the Chern-Simons term. This 
phenomenon should be observable in inelastic x- 
ray or neutron scattering experiments in solid state 
realizations or by laser assisted tunneling experi¬ 
ments on ultracold fermions in optical lattices. In 
terms of transport, the Hall conductivity clearly in¬ 
dicates the well-known plateaux ±e^/2/i for small 
electron-phonon interaction. For an increasing in¬ 
teraction, though, the plateaux are shifted away 
from each other. This is a sign that the phonons 
increase the Hall conductivity. 

A natural question is to ask whether or not 


the physics we discussed above is related to the 
fractional quantum Hall effect (FQHE), where the 
Chern-Simons field theory gives a simple and el¬ 
egant explanation of the phenomenon (Tslls^l^ . 
Within the Ginzburg-Landau approach to FQHE, 
the Chern-Simons term 


5csH = ^ (37) 

of the internal statistical electromagnetic field a 
with the flux quantum = 27r in chosen units 
and filling factor v appears in the action in order 
to impose the constrain which links a to the parti¬ 
cle density. The bosonic Ginzburg-Landau action 
should be equivalent to a fermionic action, since 
it is assumed that the quasiparticles are fermions. 
This requires a statistical field with an odd denom¬ 
inator for the corresponding filling factor: 


V = 


I 

2fc-b 1’ 


(38) 


i.e. with integer k. Including Eq. dSZl) into 
the FQHE effective action ultimately leads to the 
quantized Hall conductivity [ssll^ 


an = 1^—, (39) 

where the SI units are restored. In Eq. (l37l) . the 
Chern-Simons term appears in a general gauge. It 
is easy to see that the particular gauge choice Og = 
0 brings it to the form of Eq. da, with a different 
prefactor, though. Clearly, Eqs. (E51) and (15^ are 
analogous upon interchanging v and s/2. This fact 
suggests that, at least in the large mass limit, the 
phonon field in our model plays a role similar to the 
statistical vector potential a of FQHE. However, 
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the finite band gap of the effective action leads to tion 
the corrections of the Hall conductivity given in 
Eq. (EH). 


= -id^n, (A2) 


Appendix A: Classical equation of motion of 
the phonon field A 

Varying action ([S]) with respect to the phonon 
fields A we acquire 

where the conserved fermion current density = 
■07^V' obeys the imaginary time continuity condi- 


with n = 07o0 representing the local electron den¬ 
sity. Combination of Eq. ini and Eq. (IA2|) gives 
the classical equation of motion (constrain condi¬ 
tion) for the phonon field A coupled to the Dirac 
fermions 

^^,A^ + d^ip - 0, (A3) 

with (fi = gnj'/2, which approximately has the 
shape of the Lorenz gauge condition. In this paper 
we concentrate on quantum fluctuations around 
this classical saddle point for which the condition 
(IA3p is not valid. 


I 

Appendix B: Evaluation of the phonon loop 


Below we perform the gradient expansion in the action Eq. (ED which leads to Eqs. ED, ED and 
ED- The second term from Eq. ED reads: 


itr[^-fm] ^A^j^[0 + m] = J J G{Q)-f^G{Q + (Bl) 

with G{Q) = [—i^ + We evaluate the loop function by the method of Feynman parameter, using 


= dx 


AB Jq [(1 - x)A + xBf ■ 


(B2) 


This gives 


f fPQ [iQ+ + f) 1 


-Tr 


= -Tr / dx 


d^Q [i(l^+ m]"f[i{Q + P) 


(27r)3 ^2j 4 "Vo """7 (27r)3 [(1 - x)Q^ + x(Q + PP + 

Shifting Q ^ Q — xP symmetrizes the denominator and all terms containing odd powers of Q in the 
numerator vanish due to angular integration. The remaining expression is 

d^Q -imp -1 + Pplnlv - QluQlv + x{l - x)p-i p-f^ 


^Tr [ dx 
4 ./n 


(B3) 


(27r)3 + m'^ + x{l — x)P"^P 

Next we perform the trace. The Chern-Simons term appears from 

-tTOTr{|^7^7^} = -zTOPaTr(7„7^7^) = 2mpQe^^, 
since /r, iz = 1, 2. Next we have 

-Tr((&7^^7,.) = -9oTr(7o7;.7o7j.) - ^ 

d(l) 1 2 e d(t) 

where we made use of the angular averaging / 7 ^ QaQjS — 7:Q that 

Jo 2 Jo 

7a7^7a7,. = 7^7,.= 0. The last term simplifies to 

a^l,2 

Trp-i^p-/^ = -2po<5^i. +PaP/3Tr(7a7,x7/37,.) = -2po<^^,. + PaP/3Tr((5„^ -f 

= -2Po'^^i. + (B4) 












Hence, prior to the momentum integration we have 


dx 


d^Q {m? + - a:(l - x){P^S^^ - 2p^p^) 


(27r)3 


_|_ 77^2 _|_ _ x)P'^Y 


(B5) 


and the gradient expansion can be easily performed. To the second order in momenta it reads 


1 f d^Q + + (ll)P^ , iP^d^^-2p^p^)' 


( 27 r )2 I ((52 _|_ 77^2^2 ^Q2_^jj^2y. 


A 

— dx a;(l — x) - 


[Q^ 


,213 


[Q^ 




A§ 


where C = 1 — 


167r y/Ag + nP 
A 

1 + ( ^ 

m 


+ sgn{m)e Cpo - {P‘^5 - p p^) 


487r|r 


-O 


A2 , 

/ J 


(B6) 


which gives the loop contributions to the phonon mass Eq. m and 


both momentum dependent terms Eq. (m and (USD. For Aq 771^ we get to Eqs. (IlllllSE 


Appendix C: Kubo formula in phonon representation 


Effective actions Eq. d^D and Eq. (0 are equivalent to 

1 2 


1 


S[ip,il^,A] = -^ J d^x 


d^x + m]i; + ■ 

The fermion average in the Kubo formula Eq. (1241) can be rewritten as 


Af, + -^(V'77x^) 


(Cl) 


+\ + (^-A^ +A^ + } 

= -^(^m,^A^,o + ^ 

+ + ^=^.0 (a^ + ). 

Keeping track on Eq. dclD and assuming A^ + to bo an independent integration variable we 

realize, that the mixed terms with both phonons and fermions vanish for x ^ 0, since action dclD is 
diagonal in space. Then, the functional fermion integration in dclD can be performed which yields for 
the only remaining part Eq. (EID. _ 
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